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Chiral effective field theory can provide valuable insight into the chiral physics of hadrons when 
used in conjunction with non-perturbative schemes such as lattice quantum chromodynamics (QCD). 
In this discourse, the attention is focused on extrapolating the mass of the p meson to the physical 
pion mass in quenched QCD (QQCD). With the absence of a known experimental value, this serves 
to demonstrate the ability of the extrapolation scheme to make predictions without prior bias. 
By using extended effective field theory developed previously, an extrapolation is performed using 
quenched lattice QCD data that extends outside the chiral power-counting regime (PCR). The 
method involves an analysis of the renormalization fiow curves of the low-energy coefficients in a 
finite-range regularized effective field theory. The analysis identifies an optimal regularization scale, 
which is embedded in the lattice QCD data themselves. This optimal scale is the value of the 
regularization scale at which the renormalization of the low-energy coefficients is approximately 
independent of the range of quark masses considered. By using recent precision, quenched lattice 
results, the extrapolation is tested directly by truncating the analysis to a set of points above 380 
MeV, while temporarily disregarding the simulation results closer to the chiral regime. This tests 
the ability of the method to make predictions of the simulation results, without phenomenologically 
motivated bias. The result is a successful extrapolation to the chiral regime. 

PACS numbers: 12.39.Fc ll.lO.Jj 12.38.Aw 12.38.Gc 



I. INTRODUCTION 

In lattice quantum chromodynamics (QCD), the cal- 
culation of observables with light dynamical quarks is 
computationally intensive, and only in recent times have 
there been successful attempts to perform calculations of 
any observable at the physical point (m^ = 140 MeV) /- 
i]. Usually, some extrapolation scheme is needed if one is 
to compare theoretical calculations with the correspond- 
ing physical observables. Utilizing lattice QCD results 
spread over a larger range of quark masses naturally en- 
ables greater statistical precision in the extrapolation. 

Quenched QCD (QQCD) was introduced as a way to 
ameliorate the computational difficulty of simulating dy- 
namical fermions on the lattice. Quenched simulations 
typically have been superseded by the wide availability 
of dynamical configurations. Nevertheless, it can still be 
used as an efficient testing ground. This is particularly 
true of the chiral extrapolation problem, where the ex- 
perimentally known values may introduce a prejudice on 
a chosen form. In QQCD, the physical target point does 
not exist, and an extrapolation of moderate- mass points 
to the chiral regime provides an unbiased test of the pro- 
cedure. 

In order to discuss the chiral behaviour of the p meson 
in QQCD, one first constructs an effective field theory 
describing the relevant low-energy degrees of freedom. 
The mass of the p meson is described by a chiral ex- 
pansion in the quark mass {rriq), which includes analytic 



terms that are polynomial in to^, and non-analytic terms 
arising from chiral loop integrals. These loop integrals 
are commonly divergent, and thus it is necessary to in- 
troduce a regularization procedure. Finite-range regu- 
larization (FRR) is selected as a regularization scheme, 
which introduces a momentum cutoff scale A into the 
loop integrals. The properties of FRR allow it to be used 
with data extending outside the power-counting regime 
(PCR), at the expense of complete scheme- independence. 
As has been demonstrated, an optimal choice of regular- 
ization scale, Agcaiej can be extracted from the lattice 
simulation results [j] . A systematic uncertainty in Agcaio 
can also be estimated, which provides a range of suitable 
values for the scale obtained from the data [5| . Thus the 
scheme-dependence in using data extending outside the 
PCR can be quantified in an unbiased fashion. 



II. EXTENDED EFFECTIVE FIELD THEORY 

In chiral effective field theory (xEFT), the diagram- 
matic formulation can be used to identify the major 
contributions to the p meson mass in QQCD 0]. 
The leading-order diagrams are the double and single 
■q' hairpin diagrams as shown in Figures [T] and ^ re- 
spectively. The constant coefficients of these loop inte- 
grals are endowed with an uncertainty to encompass the 
possible effects of smaller contributions to order 0{m'^). 
Interactions with the flavour-singlet rj' are the most im- 



FIG. 1. Double hairpin tj' diagram. 
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FIG. 2. Single hairpin Vj diagram. 
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FIG. 6. Diagrammatic representation of r( propagator terms. 



finite mass — wliicli survives in the chiral limit — by 
re-summing the chain of vacuum insertions as depicted 
in Figure [Bl As a "heavy" degree of freedom, the r/ can 
then be integrated out of the of the effective field theory. 
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FIG. 3. Double hairpin quark flow diagram. 
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FIG. 4. Alternative double hairpin quark flow diagram. 
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FIG. 5. Single hairpin quark flow diagram. 



portant contributions to the p meson mass in QQCD. 
This is an artifact of the quenched approximation, where 
the r\' also behaves as a pseudo-Goldstone boson, having 
a "mass" that is degenerate with the pion. The dressing 
of the p meson by the 77' field is illustrated in Figures 
|3] through [5l Since the hairpin vertex must be a flavor- 
singlet, the mesons that can contribute are the 77' meson, 
and the w meson. The contributions from the w meson 
are insignificant due to OZI suppression and the small 
p-uj mass splitting. However, in QQCD, the 77' loop be- 
haves much as a pion loop, yet with a slightly modified 
propagator. 

In full QCD however, the 77' does not play any role 
in the low-energy dynamics. The physical 77' acquires a 



A. Loop integrals and definitions 

Using the Gell-Mann— Oakes— Renner Relation con- 
necting quark and pion masses (assuming negligible 
anomalous scaling), 777,, oc m}^ Q, the p meson mass 
extrapolation formula in QQCD can be expressed in a 
form that contains an analytic polynomial in rr?^ plus 
the chiral loop integrals (S*^): 

777,p_Q = ao + a2m\ + a477i^ 

+ £^,,,(777^, A) + S^,(7n2, A) + 0{ml). (1) 

The coefficients ai are the 'residual series' coefficients, 
which correspond to direct quark-mass insertions in 
the underlying Lagrangian of chiral perturbation theory. 
However, the non-analytic behaviour of the expansion 
arises from the chiral loop integrals. Upon renormaliza- 
tion of the divergent loop integrals, these will correspond 
with low-energy constants of the quenched xEFT. The 
extraction of these parameters from lattice QCD results 
will now be demonstrated. 

By convention, the non-analytic terms from the dou- 
ble and single hairpin integrals are Xi'^n- s-^d Xs"^^; 
spectively. The coefficients x\ and X3 of the leading- 
order non-analytic terms are scheme-independent con- 
stants that can be estimated from phenomenology. The 
low-order expansion of the loop contributions takes the 
following form: 

+ 0{ml), (2) 
S3 = &3' + btml + xtml + rat + 0{ml), (3) 



The coefficient X3 is obtained by adding the contribu- 
tions from both integrals, X3 — xi'^ +X3 ■ Each integral 
has a solution in the form of a polynomial expansion ana- 
lytic in 771^ plus non-analytic terms, of which the leading- 
order term is of greatest interest. The coefficients hi are 
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scale-dependent and therefore scheme-dependent. In or- 
der to achieve an extrapolation based on an optimal FRR 
scale, first the scale-dependence of the low-energy ex- 
pansion must be removed through renormalization. The 
renormalization program of FRR combines the scheme- 
dependent bi coefficients from the chiral loops with the 
scheme-dependent ai coefhcients from the residual series 
at each chiral order i. The result is a scheme- independent 
coefficient c,-: 



Co = flo + K"^ + bo , 

C2 = 02 + 62 + ^2 ; 

C4 = 04 + 64'' + bj , etc. 



(4) 
(5) 
(6) 



That is, the underlying coefficients undergo a renor- 
malization from the chiral loop integrals. The renormal- 
ized coefficients Ci are an important part of the extrap- 
olation technique. A stable and robust determination of 
these parameters forms the core of determining an opti- 
mal scale A'''''''''. 

The loop integrals can be expressed in a convenient 
form by taking the non-relativistic limit and performing 
the pole integration for kg. Renormalization is achieved 
by subtracting the relevant terms in the Taylor expan- 
sion of the loop integrals and absorbing them into the 
corresponding low-energy coefficients, cf. 
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(7) 
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The tilde {') denotes that the integrals are written out 
in renormalized form to chiral order 0{ml). The coeffi- 
cients Xv'v' ^^'^ Xv' related to the coefficients of the 
leading-order non-analytic terms by: 
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(9) 
(10) 



These couplings are discussed in detail below. The func- 
tion u(fc; A) is a finite-range regulator with cutoff scale 
A, which must be normalized to 1 at A:^ = 0, and must 
approach sufficiently fast to ensure convergence of the 
loop. Different functional forms of u{k; A) are equivalent 
within the PGR jl, IQ^. Different choices of u{k;A) for 
this investigation are discussed in Sec. IIIBI 

With the loop integrals specified, Eq. ^ can be rewrit- 
ten in terms of the renormalized coefficients cf. 

™P,Q = Co + C2ml + ami + S,'^,,^, (m^; A) 

+ t^,{ml;A) + 0{ml) (11) 



Co + Xi^TT + C2ml + xsm-l + c^ml 



0{ml). 



(12) 



Eq. (|lll) will be used as the extrapolation formula for 
TOp^Q at infinite lattice volume. The fit coefficients are 
Co, C2 and C4, and mp^Q is obtained by taking the square 
root of Eqs. dTTI) and p^ . It is important to note that 
the formula in Eq. is equivalent to Eq. pT|) only as 
A is taken to infinity. 

Since lattice simulations are necessarily carried out 
on a discrete spacetime, any extrapolations performed 
should take into account finite-volume effects. The low- 
energy effective field theory is ideally suited for charac- 
terising the leading infrared effects associated with the 
finite volume. In order to achieve this, each of the three- 
dimensional integrals can be transformed to its form on 
the lattice using a finite sum of discretized momenta, fol- 
lowing Armour et al. fzl], for instance: 



(2,) 



LxLyL; 



(13) 



Each momentum component is quantized in units of 
2t:/L, that is ki — ni2TT/L for integers n^. Finite- volume 
corrections (5^^^ can be written simply as the difference 
between the finite sum and the corresponding integral. 
It is known that the finite- volume corrections saturate to 
a fixed result for large values of the regularization scale 
Following the example set by this article, the value 
A' — 2.0 GeV is chosen to evaluate all finite- volume cor- 
rections independent of the FRR cutoff scale A in Eqs.Q 
and (III) . The finite- volume version of Eq. (fTTj) can thus 
be expressed: 



= Co + C2TO^ 
+ JFVC,_2 

Oiml) 



-aml + {t'^,^,{ml:A) 

cFVC/ 



{ml-A' j) + (S;;f (m^; A) + <57^(m2; A')) 



(14) 



The convention used for defining the values of xii X3j 
and the various coupling constants that occur in each, 
follows Booth For the possible different values that 
coupling constants can take, definitions by Chow & Rey 
6], Armour et al. and Sharpe [l2| are used. The 
types of vertices available are shown in Figure [71 where 
(72 and g4 occur explicitly in the two diagrams considered 
here. Booth suggests naturalness for g2 ^ 1, and that 
(74 l/Nc- These quenched coupling constants can be 
connected with the experimental value of 5tjp7r as per 
Lublinsky 13] by the relation: 



_ 1 , 
52 — 25i^P^'''^' 



(15) 



where g^^p^^ = 14 ± 2 GeV~^ and the pion decay con- 
stant takes the value /tt = 0.0924 GeV. Thus 52 is cho- 
sen to be 0.65 ± 0.09 GeV and (74 is chosen to be ap- 
proximately (72/3. The coupling between the separate 
legs of the double hairpin diagram are approximated by 
the massive constant Mq oc to^, . The next-order correc- 
tion to Mq in momentum k defines the coupling to be 
— Afg -|- Aok'^. These constants can be connected to the 
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full QCD 7]' meson mass m^/ by considering the geometric 
series of terms as previously illustrated in Figure [6] For 
the value of Mq, Booth suggest Mq 400 MeV by com- 
paring the estimate from a hairpin insertion to the result 
from the Witten-Veneziano formula 11]. In a paper by 
Duncan et al. a value of Mq ~ 900 MeV is obtained if the 
coupling constant Aq is natural. Furthermore, an analy- 
sis of the topological susceptibility leads to an estimate 
Mf) = 1.1 ± 0.2 GeV ^1^1 . In this analysis, an average 
value Mo = (400 -f 900 + 1100)/3 = 800 MeV is sensible 
as a first approximation. As a further check, consider the 
formula from Ref. [l^ |. using our normalization for the 
pion decay constant (/^ = 2 fl j^^^^^^): 



FIG. 7. Coupling types following convention introduced by Booth 
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(16) 



This formula relates the couplings Aq and Mq to 
the anomalous scaling parameter of the pion mass in 
quenched QCD, defined by: 

ml^mt'. (17) 

The parameter 5 is found to be small (and the Gell- 
Mann— Oakes— Renner Relation a good approximation), 
with a maximum value estimated by Duncan to be 
<5max — 0.03 [14]. Booth comments that the parameter 
Aq is small, and vanishes in the limit ^ oo. Nev- 
ertheless, Sharpe uses a finite value Aq ~ 0.2 [l^]. By 
using these finite values for 5 and ^o, Eq- (|T6)) leads to 
a value of Mg « 0.6 GeV^. As a result, Mq is taken to 
be 0.6 ± 0.2 GeV^ and Ao is taken to be ± 0.2. 

The coefficients Xn'ri' and Xn' can be specified in terms 
of the relevant coupling constants: 



Xri'ri' 
Xri' 



-2 m 

o 

-2 m 



929 A 



(18) 



where the couplings are defined relative to nip represent- 
ing the p meson mass in the chiral limit, which is taken 
to be 770 MeV. 



B. Finite-range regularization 

In FRR, regulator functions u{k; A) with characteristic 
scale A are inserted into the loop integrals to control the 
ultraviolet divergences that occur in the loop integrals en- 
countered. For some choices of regulator, extra regulator- 
dependent non-analytic terms arise in the chiral expan- 
sion of Eq. ([T^ . Since the correct non-analytic terms of 
the chiral expansion are regularization scale-independent 
terms, the extra non- analytic terms within working chi- 
ral order must be removed. All scale-dependence should 
be absorbed into the analytic fit parameters a^. For ex- 
ample, if a dipole regulator is chosen, the extra terms 



h^^^m%, (65^^ -I- b^^^)m% and higher-order terms occur- 
ring at odd powers of feature in Eq. p^ . One can 
avoid this by choosing a regulator that does not gener- 
ate these extra terms, up to working-order 0(rn*). Since 
the step function u^(k; A) — 9{A — k) introduces inconve- 
nient finite-volume artifacts, a 'triple-dipole' form factor 
will be chosen, defined by: 



u{k] A) 



A2 



(19) 



III. LATTICE SIMULATION DETAILS 

The calculation is performed on a 20^ x 32 lattice 
with 197 gau ge c onfigurations generated with the Iwasaki 
gauge action [1^ at /? = 2.264, and the quark propagators 
are calculated with overlap fermions and a wall-source 
technique. The lattice spacing is 0.153 fm, as determined 
from the Sommer scale parameter. 

The massive overlap Dirac operator is defined [l6l| in 
the following way so that at tree-level there is no mass 
or wavefunction renormalization 



lH: 



(20) 



where e{H) is the matrix sign function of an Hermitian 

operator H. e{H) = Hw/\Hw\ = HwKHl^Hwf'^ is 
chosen, where H-[y{x,y) — j^Dy\/{x,y). Here D^/ is the 
usual Wilson-Dirac operator on the lattice, except with 
a negative mass parameter —p = I /2k — 4 in which Kc < 
K < 0.25. Ta king k — 0.19 in the calculation corresponds 
to p = 1.368 [11 [H. 

In Figure[T2the simulation results for the vector meson 
mass are shown for a range of quark masses. 

The data displayed in Figure[T2]are split into two parts. 
All the data left of the solid vertical line is unused for ex- 
trapolation and kept in reserve. Indeed, the authors per- 
forming the extrapolation were blind to these data. This 
is so that the extrapolation can be checked against these 
known data points once the extrapolation is established. 
In other words, the results of the chiral extrapolation are 
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TABLE I. Quenched lattice QCD data for the p meson mass rup at 
various pion mass squared values m^. The statistical uncertainty 
of the is negligible. The values of the bare quark mass rUq are 
also included for comparison. The lattice size is 20^ X 32, with 
a lattice spacing of 0.153 fm. Entries below the line (underneath 
= 0.143 GeV^) remained hidden until the extrapolation was 
determined. 
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FIG. 8. (color online). Effective mass plot corresponding to the 
simulation at = 0.143 GeV^. Only the wall source results are 
plotted. The point source results are not used in the analysis. 
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FIG. 9. (color online). Effective mass plot corresponding to the 
simulation at = 0.080 GeV'^. Only the wall source results are 
plotted. The point source results are not used in the analysis. 
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genuine predictions of the hidden lattice results. Only 
the data points to the right of the solid vertical line are 
used for extrapolation. The full set of data is also listed in 
Table m which also includes the bare quark mass values. 
In addition, effective mass plots corresponding to four 
lighter pion masses are included, in Figures [8] through 

EI 

To estimate finite-volume effects using overlap 
fermions, quenched lattices of volumes 16'^ x 28 and 
X 28 with a — 0.2 fm are used. For a pion mass 
of 180 MeV, mpsL w 3, and the finite- volume correction 
is approximately 2.7 MeV: about 1.5% of the pion mass 
pl |. The current 20^ x 32 lattice with a = 0.153 fm is 
about the same physical size as that of a 16^ x 28 lattice 
and a similar finite- volume correction is expected. To es- 
timate the finite- volume correction of the lowest p meson 
mass at « 200 MeV, the same percentage of error 
is used, and a shift of SltUp « 13 VIeV to the p mass 



FIG. 10. (color online). Effective mass plot corresponding to the 
simulation at = 0.053 GeV^. Only the wall source results are 
plotted. The point source results arc not used in the analysis. 



FIG. 11. (color online). Effective mass plot corresponding to the 
simulation at = 0.041 GeV^. Only the wall source results are 
plotted. The point source results are not used in the analysis. 



6 



is calculated for the p meson mass of rup w 917 MeV. 
This is about half of the statistical error of the lattice 
data. It should be noted that the data that will be used 
in chiral extrapolations are those with pion mass greater 
than 400 MeV, with rapsL > 6.2. The predictions are 
extended to the region with pion mass less than 400 MeV 
and compared with the lattice data. 

With regards to possible lattice artifacts, the lattice 
results analyzed are based on the overlap fermion on 
quenched gauge configurations at one lattice spacing. 
Even though the overlap fermion has relatively smaller 
0{a) errors, the 0(0^) correction toward the continuum 
limit has not been taken into account. With a spatial 
size of 3.06 fm, mTfO ~ 3 for the smallest pion mass at 
200 MeV is somewhat smaller than mTrO = 4, beyond 
which the finite volume effect has been considered to be 
small. For mTrO ~ 3, the previous study described in 
Ref. [3l estimates that the finite-volume correction is 
approximately 3% which is smaller than the statistical 
error of the pion mass. 

The enhancement of zero modes effects in QQCD pri- 
marily affects the pseodoscalar and scalar mesons. Since 
all the zero modes appear in one chiral sector in each 
gauge configuration, the pseodoscalar and scalar mesons 
will have a leading singularity from the zero modes. 
These appear in both the quark and antiquark propaga- 
tors in the meson correlator Nevertheless, the vec- 
tor and axial vector mesons have only a l/m singularity, 
which is a less dramatic effect. In either case, the quan- 
tity that determines the size of the zero mode effects is 
m'EV in the p-regime (20j . It has been demonstrated 
that when mY,V > 5, the zero mode effect is hardly de- 
tectable [isl . [2lj | . For all pion masses displayed in Figure 
[T^ mTiV > 7. Therefore, there is no reason to suggest 
that there is a zero mode contribution to the p meson 
correlators being studied. 



IV. EXTRAPOLATION RESULTS 
A. Renormalization flow curves 

In order to produce an extrapolation to each test value 
of mj, a finite-range regularization scale A must be se- 
lected. As an example, one can choose a triple-dipole 
regulator at A = 1.0 GeV. By using Ea. (IT4l) . finite- and 
infinite-volume extrapolations are shown in Figure 1131 
Note that the values selected for the finite- volume 
extrapolations exactly correspond to the 'missing' low- 
energy data points set aside earlier. The physical point 
= 0.0196 GeV^ is included as well. 

Now the regularization scale-dependence of low-energy 
coefficients cq, C2 and C4 is investigated for various upper 
limits of the range of pion masses. The renormalization of 
these low-energy coefficients is considered for a series of A 
values. The aim is to obtain renormalization fiow curves, 
each corresponding to a different value of maximum pion 
mass, m?. Thus the behaviour of the renormalization 
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FIG. 12. (color online). Quenched lattice QCD data for tlie p 
meson mass. The dashed vertical line indicates the physical pion 
mass and the solid vertical line shoves how the data set is split into 
two parts. The lower-mass portion of the data was not known at 
the time of extrapolation. 
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FIG. 13. (color online). A test extrapolation based on the 
four lightest original data points (excluding the low-energy set) 
as shown. Both the finite- and infinite-volume results are shown 
for a triple-dipole regulator at A = 1.0 GeV. The dashed vertical 
line indicates the physical pion mass. 



of the low-energy coefficients can be examined as lattice 
data extend further outside the PGR. Figures [T4l through 
1161 show the renormalization flow curves for each of cq, 
C2 and C4. Note that each data point plotted has an 
associated error bar, but for the sake of clarity only a 
few points are selected to indicate the general size of the 
statistical error bars. Using the procedure described in 
Ref. [4], the optimal regularization scale is identified by 
the value of A that minimizes the discrepancies among 
the renormalization flow curves. This indicates the value 
of regularization scale at which the renormalization of 
Co, C2 and C4 is least sensitive to the truncation of the 
data. Physically, this value of A can be associated with 
an intrinsic scale related to the size of the source of the 
pion cloud. 

By examining Figures [14] through I16[ increasing 
"^^.max leads to greater scheme-dependence in the renor- 
malization, since the data sample lies further from the 
PGR. Gomplete scheme- independence would be indicated 
by a horizontal line at the physical point. Since the ef- 
fective field theory is calculated to a finite chiral order, 
complete scheme-independence across all possible values 
of A will not occur in practice. Note that an asymptotic 
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FIG. 14. (color online). Behaviour of cq vs. A. A few points are 
selected to indicate the general size of the statistical error bars. 
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FIG. 15. (color online). Behaviour of C2 vs. A. A few points are 
selected to indicate the general size of the statistical error bars. 
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FIG. 16. (color online). Behaviour of C4 vs. A. A few points are 
selected to indicate the general size of the statistical error bars. 



value is usually observed in the renormalization flow as A 
becomes large, indicating that the higher-order terms of 
the chiral expansion are effectively zero. However, these 
asymptotic values of the low-energy coefficients are poor 
estimates of their correct values, as previously demon- 
strated in a pseudodata model Instead, the best 
estimates of the low-energy coefficients lie in the iden- 
tification of the intersection point of the renormalization 
flow of the low-energy coefficients. It is also of note that, 
for small values of A, FRR schemes break down. The reg- 
ularization scale must be at least large enough to include 
the chiral physics being studied. 



B. Optimal regularization scale 

The optimal regularization scale A'^'^'^''' can be obtained 
from the renormalization flow curves using a chi-square 
analysis described below. In addition, the analysis will 
allow the extraction of a range for A*"^^'°. Knowing how 
the data are correlated, the systematic uncertainties from 
the coupling constants and A'^'^^''' will be combined to ob- 
tain an error bar for each extrapolation point. Of partic- 
ular interest are the values of rup^q at the values of 
explored in the lattice simulations but excluded in the 
chiral extrapolation. 

To obtain a measure of the uncertainty associated with 
an optimal regularization scale, a Xdof function is con- 
structed. This function should allow easy identification of 
the intersection points in the renormalization flow curves, 
and a range associated with this central regularization 
scale. The flrst step is to plot Xdof against a series of A 
values. The relevant data are the extracted low-energy 
coefficients with differing values of „jax- ^ plot of 
Xdof i^ constructed separately for each renormalized co- 
efficient c (with uncertainty 6c): 



\dof 



_^^(c(z^A^^c^W 
(<5c(*;A))2 



1 ^ 

1=1 



(21) 



2 

-7r,niax 



for i corresponding to fits with differing values of m 
(n = 8). The theoretical value C"^ is given by the weighted 
mean: 



ELic(z;A)/(fc(»;A)f 

E;=ii/(Mj;A)f 



(22) 



The Xdof plots using a triple-dipole regulator are shown 
in Figures [T7] through 1191 The optimal regularization 
scale A'^'^'^^'' is taken to be the central value Acontrai of 
each plot. The upper and lower bounds obey the condi- 
tion xlof < xlof,min + ^/(dof). The results for the opti- 
mal regularization scale and the upper and lower bounds 
are shown in Table |lll It is remarkable that each low- 
energy coefficient leads to the same optimal value of A, 
i.e. Accntrai = 0.67 GeV. By averaging the results among 
Co, C2, and C4, the optimal regularization scale Agcaic for 
the quenched p meson mass can be calculated for this 



data set: Agcaio = 0.67 



+0.09 
-0.08 



GeV. 



The result of the final extrapolation, using the estimate 
of the optimal regularization scale Agcaic — 0.67lg'Qg 
GeV, and using the initial data set to predict the low- 
energy data points, is shown in Figure The extrap- 
olation to the physical point obtained for this quenched 
data set is: ^^^^(to^ pj^^J = 0.925^;^:^^^ GeV, an uncer- 
tainty of less than 6%. 

Note that each extrapolation point displays two error 
bars. The inner error bar corresponds to the systematic 
uncertainty in the parameters only, and the outer error 
bar corresponds to the systematic and statistical uncer- 
tainties of each point added in quadrature. Also, the 
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FIG. 17. (color online) 
the renormalization flow 
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FIG. 18. (color online). Xdof "^2 versus A, corresponding to 
the renormalization flow curves displayed in Figure [TS] 
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the renormalization flow 
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TABLE II. Values of the central, upper and lower regularization 
scales, in GeV, obtained from the Xdof ™^lysis of cq, C2 and C4, 
displayed in Figures 1 171 through 1191 



scale (GeV) 


CO (FigdZl) 


C2 (FigdU) 


C4 (FigdU) 


-^central 


0.67 


0.67 


0.67 


-^uppcr 


0.78 


0.75 


0.75 


-Slower 


0.58 


0.59 


0.60 



infinite- volume extrapolation curve is displayed in order 
to illustrate the effect of finite- volume corrections to the 
loop integrals. 

In Figure 1211 the extrapolation predictions are com- 
pared against the actual simulation results, which were 
not included in the fit. Note that both the extrapolations 
and the simulation results display the same non-analytic 
curvature near the physical point. Figure [22] shows the 
data plotted with error bars correlated relative to the 
lightest data point in the original set, = 0.143 GeV^. 
This aids in clarifying the plot from Figure [5l] by re- 
moving much of the correlated statistical error in the 
lattice data, and allows us to be even more stringent in 
determining whether the extrapolation is successful. It 
is notable that the extrapolated results are consistent 
with the lattice data even after having removed the cor- 
related statistical error. To highlight the importance of 
this application of an extended xEFT, a simple linear 
fit is included in Figure [221 By ignoring low-energy chi- 
ral physics, the linear fit is statistically incorrect at the 
physical point. Note also that all of the missing original 
data points are consistent within the extrapolations' sys- 
tematic uncertainties. After statistical correlations are 
subtracted, the extrapolated points correspond to an er- 
ror bar almost half the size of that of the lattice data 
points. In order to match this precision at low energies, 
the time required in lattice simulations would increase by 
approximately four times. 

In order to check if scheme-independence is recovered 
using data within the PGR, the low-energy data that 
were initially excluded from analysis can now be treated 
in the same way. That is, renormalization flow curves can 
be constructed as a function of A for sequentially increas- 
ing The results are shown in Figures [23l through 
[211 Clearly, the renormalization fiow curves for each plot 
corresponding to cq, C2 and C4 are flatter than those of 
the initial analysis, indicating a reduction in the regular- 
ization scale-dependence due to the use of data closer to 
the PGR. One is not able to extract an optimal regular- 
ization scale from these plots, as shown in the behaviour 
of Xdo/' displayed in Figures [HI through [2H1 However, 
each Xdof curve provides a lower bound for the regular- 
ization scale, where FRR breaks down [3], as discussed 
in Section HVAI These lower bounds are: A^"^^^ = 0.39 
GeV, Ag^^^ = 0.52 GeV and A^^^^ = 0.59 GeV. 

The statistical error bars of the low-energy coefficients 
corresponding to a small number of data points in Fig- 
ures [221 through [221 is large, and a statistical difference 
among the curves does not appear until m-^ max ~ 0.11 
GeV^. Thus the identification of an optimal regulariza- 
tion scale will be aided by incorporating data correspond- 
ing to even larger values of to J ^ax- By considering all 
of the available data, the behaviour of Xdof^ displayed 
in Figures [^ through 13 1 [ resolve precise optimal regular- 

0.71 GeV 



ization scales: A;;^^^^^! 



0.72 GeV, A^=^,,,j 



and A^*jjjj.g^j = 0.71 GeV. The systematic errors obtained 
from each Xdof curve seem arbitrarily constrained as a 
consequence of including more data points, which extend 
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FIG. 20. (color online). Extrapolation at As^alc = CeTlpiog '^'^^ 
based on Kentucky Group data, and using the optimal number of 
data points, corresponding to m^ max = 0.35 GeV^. The inner 
error bar on the extrapolation points represents purely the sys- 
tematic error from parameters. The outer error bar represents the 
systematic and statistical error estimates added in quadrature. 
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FIG. 21. (color online). Comparison of chiral extrapolation predic- 
tions (blue diamond) with Kentucky Group data (red cross). Ex- 
trapolation is performed at Agcalo = 0-67^Q g| GeV, and using the 
optimal number of data points, corresponding to ^^-^ = 0.35 
GeV^. The inner error bar on the extrapolation points represents 
purely the systematic error from parameters. The outer error bar 
represents the systematic and statistical error estimates added in 
quadrature. 
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FIG. 22. (color online). Comparison of chiral extrapolation pre- 
dictions (blue diamond) with Kentucky Group data (red cross), 
with errors correlated relative to the point at = 0.143 GeV^. 
This is done simply to clarify the plot in Figure 1211 by removing 
much of the correlated statistical error. Extrapolation is performed 



original data 
fin.vol. extrap 
naive linear fit 



at A=, 



: 0.67"^ 



GeV, and using the optimal number of data 



points, corresponding to "i^ max = 0.35 GeV^. The error bar on 
the extrapolation points represents the systematic error only. A 
simple linear fit, on the optimal pion mass region, is included for 
comparison. 
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FIG. 23. (color online). Behaviour of cq vs. A including the ini- 
tially excluded low-energy data. A few points are selected to indi- 
cate the general size of the statistical error bars. 
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FIG. 24. (color online). Behaviour of C2 vs. A including the ini- 
tially excluded low-energy data. A few points are selected to indi- 
cate the general size of the statistical error bars. 
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FIG. 25. (color online). Behaviour of C4 vs. A including the ini- 
tially excluded low-energy data. A few points are selected to indi- 
cate the general size of the statistical error bars. 



well outside the chiral regime, and possibly outside the 
applicable region of FRR techniques. This issue is ad- 
dressed in the ensuing section. 



C. Optimal pion mass region 



In this section, a robust method for determining an 
optimal range of pion masses is presented. This range 
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FIG. 26. (color online). Xdof '^"'^ '^O versus A, corresponding to 
the renormalization flow curves displayed in Figure 1231 A lower 
bound for the regularization scale is found: A^^^^^ = 0.39 GeV. 
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FIG. 27. (color online). Xdof^ ^'-'^ versus A, corresponding to 
the renormalization flow curves displayed in Figure 1241 A lower 



bound for the regularization scale is found: A, 



0.35 



: 0.52 GeV. 
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FIG. 28. (color online). Xdof '^"'^ versus A, corresponding to 
the renormalization flow curves displayed in Figure 1251 A lower 
bound for the regularization scale is found: Af*' = 0.59 GeV. 



corresponds to an optimal number of simulation results 
to be used for fitting. First, consider the extrapolation 
of the quenched p meson mass, which can now be com- 
pleted. The statistical uncertainties in the values of cq, 
C2, C4 are dependent on rn^ As a consequence, the 
uncertainty in the extrapolated p meson mass m°^^ must 
also be dependent on Since the estimate of the 

statistical uncertainty in an extrapolated point will tend 
to decrease as more data are included in the fit, one might 
naively choose to use the largest j^^x value possible 
in the data set. However, at some large value of ^g^^j 
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FIG. 29. (color online). Xdof '^"'^ "^0 versus A, corresponding to 
all available data, including the low-energy set. 
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FIG. 30. (color online). Xdof ^'-'^ "^2 versus A, corresponding to 
all available data, including the low-energy set. 
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FIG. 31. (color online). Xdof "^4 versus A, corresponding to 
all available data, including the low-energy set. 



FRR xEFT will not provide a valid model for obtaining 
a suitable fit. At this upper bound of applicability for 
FRR xEFT, the uncertainty in an extrapolated point is 
dominated by the systematic error in the underlying pa- 
rameters. This is due to a greater scheme-dependence 
in extrapolations using data extending outside the PGR, 
meaning that the extrapolations are more sensitive to 
changes in the parameters of the loop integrals. Thus 



there is a balance point to^ , 



^, where the sta- 



tistical and systematic uncertainties (added in quadra- 
ture) in an extrapolation are minimized. 

In order to obtain this value to J consider the be- 
haviour of the extrapolation of the p meson mass to the 



physical point '7i°'g(TO^ pj^y^), as a function of 
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FIG. 32. (color online). Behaviour of the extrapolation of the 



quenched p meson mass to the physical point m^Q{m' 
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max using the initial data set, which excludes the lowest mass 
data points. In each case, cq is obtained using the scale A^entral i^or 
a triple-dipole regulator) as obtained from the Xdof a-^alysis. The 
error bars include the statistical and systematic uncertainties in cq 



added in quadrature. The optimal value rhj. , 
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FIG. 33. (color online). Behaviour of the extrapolation of the 
quenched p meson mass to the physical point "^p^q('ti^ phys) 

max using the complete data set, which includes the lowest mass 
data points. In each case, cq is obtained using the scale A^entral (for 
a triple-dipole regulator) as obtained from the Xdof a-ualysis. The 
error bars include the statistical and systematic uncertainties in cq 
added in quadrature. The optimal value ^^-^ = 0.20 GeV^. 



Treating the parameters: A'"^'*''', g2, 54, Mq and ^0 as 
independent, their systematic uncertainties from these 
sources are added in quadrature. In addition, the sys- 
tematic uncertainty due to the choice of the regulator 
functional form is roughly estimated by comparing the 
results using the double-dipole and the step function. 
These functional forms are the two most different forms 
of the various regulators considered, since the dipole was 
excluded due to the extra non-analytic contributions it 
introduces. The results for the initial and complete data 
sets are shown in Figures EH and [331 respectively. Note 
that the systematic uncertainty due to A^'^'^^'' is included 
for chiral order O(m^). 

Figure |32] indicates an optimal value ^^^^ = 0.35 
GeV^, which will be used in the final extrapolations, in 
order to check the results of this method with the low- 
energy data. By using only the data contained in the 
optimal pion mass region, constrained by an es- 



timate of the optimal regularization scale may be cal- 
culated with a more generous corresponding systematic 
uncertainty. The value A'^'^'^'® = 0.64 GeV is the average 
of Aff'', All"^" and Aff"" using this method. The xL/ 
analysis does not provide an upper or lower bound at 
this value of mi ,„„^. Note that these two estimates of 
the optimal regularization scale are consistent with each 
other. Both shall be used and compared in the final anal- 
ysis. Figure [33l indicates an optimal value m^^^^^ = 0.20 
GeV^ for the complete data set, which includes the low- 
energy data. A higher density of data in the low-energy 
region serves to decrease the statistical error estimate of 
extrapolations to the low-energy region. The correspond- 
ing value of A'''^'^'° is unconstrained in this case, since the 
data lie close to the PGR. The breakdown of the system- 
atic error bar into its constituent uncertainties is listed 
in Table Uni 

The values of cq, C2 and C4 for both the original data 
set and the complete data set are shown in Table IIV| 
with statistical error estimate quoted first and systematic 
uncertainty due to the parameters A'^'^^^'^, 92, 54, -^0 , ^Oj 
and the regulator functional form quoted second, in this 
order. In the case of the original data set, the value of C4 
is not well-determined, due to the small number of data 
points used. In the case of the complete data set, the 
results are dominated by statistical uncertainty and also 
results in an almost unconstrained value of C4. Even if 
^scaio jg quite well-determined, as observed in Figures fTTl 
through [ini the value of C4 itself is sensitive to changes in 
the regularization scale A, as evident from Figure[T6l The 
coefficients of the complete set are less well-determined 
due to the fact that m^^max = 0-20 GeV^, leaving only 
low-energy results with large statistical uncertainties for 
fitting. 

The result using the estimate of the optimal regular- 
ization scale A'^'^'^''' = 0.64 GeV, with the systematic un- 
certainty calculated by varying A across all suitable val- 
ues, and using the initial data set, is shown in Figure l34l 
The extrapolation to the physical point obtained for this 
quenched data set is: mp^glm^ pj^yj = 0.922t°:°^|5 GeV, 
an uncertainty of approximately 7%. Figure [3S] shows 
the data plotted with error bars correlated relative to 
the hghtest data point in the original set, = 0.143 
GeV^, using Ascaie = 0.64 GeV, and varying A across its 
full range of values. This naturally increases the estimate 
of the systematic uncertainty of the extrapolations, but 
also serves to demonstrate how closely the results from 
lattice QGD and xEFT match. 



CONCLUSION 



A technique for isolating an optimal regularization 
scale, established in Ref. M!], was tested in quenched 
QCD through an examination of the quenched p meson 
mass. The result is a successful extrapolation based on 
an extended effective field theory procedure. By using 
quenched lattice QCD results that extended beyond the 
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TABLE III. Results for the quenched p meson mass for different values of maxi extrapolated to the physical point, corresponding to 
Figures 1321 (for the original data set) and 1331 (for the complete data set). The uncertainty in rn''^Q{m'^ phys) provided in the following 
order: the statistical uncertainty, the systematic uncertainty due to the intrinsic scale, g2, Qi, M^, Aq and the regulator functional form, 
respectively. 





^7^Q{ml,phyJ (GeV): original set 


"ip^Q ("T-^,phyJ (GeV): complete set 


0.059 


- 


0.956(234)(1)(0)(0)(1)(0)(0) 


0.068 


- 


0.938(139)(1)(1)(0)(1)(0)(0) 


0.080 


- 


0.934(87)(1)(1)(0)(1)(0)(0) 


0.094 


- 


0.939(64)(2)(1)(0)(1)(0)(1) 


0.114 


- 


0.930(47)(3)(2)(0)(2)(0)(0) 


0.143 


- 


0.932(34)(5)(3)(0)(4)(1)(0) 


0.172 


- 


0.929(31)(6)(4)(0)(5)(1)(0) 


0.204 


- 


0.929(29)(9)(5)(0)(6)(1)(0) 


0.241 


0.930(110)(27)(14)(0)(17)(4)(6) 


0.927(27)(12)(7)(0)(9)(2)(0) 


0.288 


0.899(62)(31)(1)(0)(1)(1) 


0.933(24)(17)(10)(0)(12)(3)(1) 


0.347 


0.922(43)(37)(11)(0)(13)(28)(17) 


0.932(23)(21)(11)(0)(13)(3)(4) 


0.422 


0.948(29)(45)(23)(1)(28)(7)(8) 


0.930(20)(29)(14)(0)(16)(4)(8) 


0.515 


0.932(23)(51)(19)(1)(23)(6)(19) 


0.929(19)(37)(15)(0)(18)(4)(13) 


0.610 


0.921(18)(63)(18)(1)(22)(5)(25) 


0.926(16)(54)(18)(0)(21)(5)(22) 


0.676 


0.915(12)(74)(18)(1)(22)(6)(32) 


0.926(14)(66)(19)(1)(23)(6)(25) 


0.743 


0.919(13)(79)(22)(1)(26)(7)(29) 


0.922(12)(74)(21)(1)(26)(7)(27) 


0.867 


0.933(9)(103)(32)(1)(39)(12)(35) 


0.923(8)(100)(29)(1)(35)(11)(38) 


1.062 


0.928(7)(115)(32)(1)(38)(12)(43) 




1.212 


0.926(7)(121)(31)(1)(37)(12)(38) 




1.329 


0.921(6)(132)(31)(1)(37)(12)(43) 




1.742 


0.915(5)(146)(30)(1)(37)(13)(48) 




2.187 


0.910(4)(175)(30)(1)(37)(14)(55) 




3.150 


0.902(3)(197)(29)(1)(36)(14)(61) 





TABLE IV. The values of cq, C2 and C4 as obtained from both the original data set and the complete set, which includes the low-energy 
data. In each case, the coefficients are evaluated using the scale Acontral (for a triple-dipole regulator) as obtained from the Xdof a-n^lysis. 
The value of m^^^^-^ used is that which yields the smallest error bar in adding statistical and systematic uncertainties in quadrature. For 
the initial data set, ^^-^ = 0.35 GeV^. For the complete data set, ^^-^ = 0.20 GeV^. The statistical uncertainty is quoted in the 
first pair of parentheses, and the systematic uncertainty is quoted due to the parameters, in the following order: A'^'^^''^, 32, 94, Mq, Aq 
and the regulator functional form. 



co(GeV2) 


C2 


C4(GeV-2) 


original set 1.31(5)(10)(4)(0)(5)(4)(8) 
complete set 1.35(4) (3) (36) (16) (60) (166) (113) 


7.9(4)(25)(2)(0)(2)(1)(4) 
6.8(5)(17)(13)(1)(17)(14)(11) 


-16.2(7)(382)(3)(0)(3)(1)(4) 
-3.3(16)(359)(23)(1)(28)(12)(1) 



power-counting regime, an optimal regularization scale 
was obtained from the renormalization flow of the low- 
energy coefhcients. 

An optimal value of the maximum pion mass to be 
used for fitting was also calculated, and this resulted in 
an alternative estimate of the value of the optimal regu- 
larization scale, which was consistent with the first result. 
The mass of the p meson was calculated in the low-energy 



region, including the physical point, using each estimate 
of the optimal regularization scale, and both results were 
compared. The results of extrapolations using xEFT, 
and the results of lattice QCD simulations, were demon- 
strated to be consistent. The extrapolation correctly pre- 
dicts the low-energy curvature that was observed when 
the low-energy lattice simulation results were revealed. 

In full QCD, using dynamical fermions, the process 
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FIG. 34. (color online). Comparison of chiral extrapolation pre- 
dictions (blue diamond) with Kentucky Group data (red cross). 
Extrapolation is performed at Alcaic = 0.64 GeV, varied across 
the whole range of A values, and using the optimal number of data 
points, corresponding to = 0.35 GeV^. The inner error bar 

on the extrapolation points represents purely the systematic error 
from parameters. The outer error bar represents the systematic 
and statistical error estimates added in quadrature. 
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FIG. 35. (color online). Comparison of chiral extrapolation pre- 
dictions (blue diamond) with Kentucky Group data (red cross), 
with errors correlated relative to the point at = 0.143 GeV^. 
This is done simply to clarify the plot in Figure 1341 by removing 
much of the correlated statistical error. Extrapolation is performed 
f't Agcale = 0.64 GeV, varied across the whole range of A values, 
and using the optimal number of data points, corresponding to 
max ~ 0-35 GeV'^. The error bar on the extrapolation points 
represents the systematic error only. A simple linear fit, on the 
optimal pion mass region, is included for comparison. 



/9 — >■ TTTT contributes to the p meson mass. This means 
that near the chiral Hmit, the tttt component of the p nec- 
essarily involves a hard momentum scale, and therefore 
is not amenable to the standard methods of low-energy 
expansions, as entailed by xPT. Therefore, one needs to 
resort to alternative techniques in such instances. 



However, since there exists no experimental value for 
the mass of a particle in the quenched approximation, 
this analysis demonstrates the ability of the technique to 
make predictions without phenomenologically motivated 
bias. The results clearly indicate a successful procedure 
for using lattice QCD data outside the power-counting 
regime to extrapolate an observable to the chiral regime. 
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